Quantum Entanglement and fixed point Hopf bifurcation 
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We present the qualitative differences in the phase transitions of the mono-mode Dicke model 
in its integrable and chaotic versions. We show that a first order phase transition occurs in the 
integrable case whereas a second order in the chaotic one. This difference is also reflected in the 
classical limit: for the integrable case the stable fixed point in phase space suffers a bifurcation of 
Hopf type whereas for the second one a pitchfork type bifurcation has been reported. 
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I. INTRODUCTION 



Entanglement is a fundamental characteristic of quan- 
tum systems which has lately provided for impressive 
progress in several areas such as quantum information P| , 
quantum cryptography Q and teleportation In this 
context, several measures of this property have been pro- 
posed, such as entropy-like quantities 0|, negativity 0, 
concurrence and so forth. In particular, one expects 
strong correlations to be present in quantum phase tran- 
sitions (QPT) 7]. Therefore, in this context it seems par- 
ticularly interesting to study the interplay between quan- 
tum correlations and the appearance of a QPT. Since the 
quantification of the degree of entanglement (or correla- 
tions) in a system is not unique, in recent years several ex- 
amples of specific systems using different measures have 
been put forth. More specifically, ground state correla- 
tions have been investigated as a function of a coup ling 
parameter. In references |llS[l3[IllilEllMC|| 
the ground state two atom concurrence has been stud- 
ied and shown to exhibit a maximum at the transition 
point, in the context of some spin-spin or spin-boson 
models . As pointed out recently by Hines-McKenzie- 
Milburn (HMM) jl^, in the classical regime correspond- 
ing to the qualitative change in the quantum ground 
state it corresponds a change in minimum energy sta- 
ble fixed point in phase space. As the parameter value of 
the Hamiltonian is changed, a quantum instability of the 
ground state is followed (at A = Ac), and a qualitative 
change in phase space structure of the system. In their 
work, the pitchfork (emergence of two new stable fixed 
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point) type bifurcation has been focused and associated 
to a spike of the entanglement in the mean field limit 
(N-^oo). 

The purpose of the present work is firstly to add to 
these previous contributions a study of an integrable ver- 
sion of Dicke's model [TsL IT^ . and to show that the 
phase transition in this case is radically different from 
the one in the non-integrable situation as revealed by the 
linear entro py ( adopted here as a measure of quantum 
correlations) poj. Reference 0| proposes, based on spe- 
cific model studies how the entanglement in a nonlinear 
bipartite system can be associated with a fixed point bi- 
furcation in the classical dynamics. This conjecture thus 
contemplates pitchfork (one-to-two) bifurcation. Our re- 
sults here suggest that first order transitions can present 
a different type of bifurcation. We show that in the inte- 
grable version of the model the abrupt change in entan- 
glement content of the ground state (GS) is associated 
to a Hopf type bifurcation (one-to- infinitely many de- 
generate fixed points). In fact this also occurs for the 
Jahn- Teller |21| and dimer models. So, we suggest 
that when the QPT is first order the corresponding clas- 
sical instability is of a Hopf type bifurcation. This idea is 
supported by the analysis of the GS Wigner function in 
the integrable and non-integrable situations. The max- 
ima of these functions exactly follow the classical fixed 
points. 



II. THE MODEL AND THE BEHAVIOR OF 
THE ENTROPY 

The quantum Hamiltonian we use is written in the 
form 

H^Ho + Hg + h'g (1) 
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with 



Ho 

Hg = 



TiLJO^a + fltJ: 

G 



G' 



aJ+ + a) J- 



(2) 



where we consider in ah numerical calculations shown 
uj — e — 1 (resonant case). When G = G' we recover the 
usual single mode Dicke Hamiltonian. The advantage 
of working with Hamiltonian is that the integrable 
regime can be obtained just by setting G' = (or G = 0), 
and it becomes easier to explore various mixed regimes. 
The Super-radiant phase transition is present in all these 
situations at G + G' = e HI. 
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Figura 1: Linear entropy S versus \ = G/e for parameter 
values J = 1.5,4.5,7.5, h = l,uj — e = 1 : integrable case 
(G' — 0), with a clear transition at Ac = 1. 



A, the smaller the steps become, and finally a plateau 
is attained for large values of the coupling. The plateau 
values of the entropy increases with increasing value of 
J, indicating a larger number of basis states participat- 
ing in the entanglement of the two subsystems. At the 
limit J >> 1 we expect the height of the steps to remain 
of the same size, but the width to become shorter and 
shorter, whereas the plateau tends to 1. In this case, the 
classical analog of this behavior stems from the fact that 
in the classical limit the equilibrium points present bifur- 
cations at the same parameter values as in the quantum 
case. The bifurcations as shown in [2^, in both - inte- 
grable and non-integrable - cases are also very distinct. 
The classical limit is discussed below. 



III. THE CLASSICAL ANALOG AND BIFURCATION 
OF EQUILIBRIA 

The classical Hamiltonian corresponding to eq. ^ can 
be obtained by a standard procedure using the coherent 
states of spin and boson 25]: 



with 



\w) = {1 + wwy' e'^'^+ \J,-J) 



\v) = e-™/2e"''+ |0) 



Pi + iqi 



^AJ-ipi + qfY 



(3) 



(4) 



(5) 



The non-inte grab le GS features of this model have been 
studied in ref. |l4| using the concurrence as a measure 
of GS correlations. It is shown that concurrence attains 
a maximum when A — A™ (N) . This maximum is a func- 
tion of the number of atoms and tends to the critical 
value Ac — Am(oo) as A'^ ^ cxd. In the present work 
we use the linear entropy of the atomic subsystem as a 
measure of the quantum correlations. We observe the 
same qualitative behavior which is therefore not shown 
here. However, for the integrable situation G' ~ shown 
in FigQ] the atomic linear entropy (ALE) is plotted as 
a function of A = G/e for N = 2J ~ 3,9,15 atoms. 
The ALE presents an abrupt jump from zero to 0.5 at 
Ac = G/e — 1 for all values of J. Moreover, the en- 
tropy has a monotonic increase in steps of finite size: for 
example, before the first jump (A < Ac) only one ba- 
sis state contributes to the ground state. It is a well 
known fact that below A = Ac = 1 the ground state of 
the integrable version of the model is an isolated basis 
state \n = 0,m ~ — J)- It is not connected to any oth- 
er basis state via the above interaction for G' = [l^ . 
until the interaction term of the Hamiltonian becomes 
dominant and the ground state becomes a state of the 
form (|1, — J) -|- |0, — J -I- 1)), which is orthogonal to 
the previous fundamental state. The larger the value of 



{P2 + iq2) 



(6) 



The classical Hamiltonian is defined as {wv \ H\ wv) |2y|, 
and in this case there results a nonlinearly coupled two 
degrees of freedom Hamiltonian function in terms of the 
parametrization (O and (|SJ): 



^(«l,Pl,92,P2) 



4J 



{G+P1P2 + G-qiq2) 



(7) 



where G± — G ± G'. The integrable situation corre- 
sponds to G' = and G± ~ G. 

The stable fixed points corresponding to the equilibri- 
um positions in phase space are defined via Hamilton's 
equations: 



91 = 
Pi 



-epi 



%(2J-i70^/^ 

V ZJ 



2J(2J-i/i)i/2 



{G+P1P2 + G_gig2) = (8) 
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+ 



Pi = -egi 
qi 



2J(2J-i7i)i/2 
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(G+P1P2 + G-qiq2) = (9) 



up,~^{2J-H,f/^^Q (10) 



2J 



P2 = t^g2 



^(2J-Hi)^/^ = (11) 
V 2 J 



where i/i = ^itE?. 



In ref. ||2Jj it has been shown that there are two solu- 
tions for: 

A. Integrable case: 

i) trivial solution (stable for G/e < 1) : 



91 Pi = 92 = P2 = 0, (the origin) 
ii) non-trivial solution (stable for G/e > 1): 



Rn 



2 I 2 
12 +P2 



J 



G* 



G2w2 



(12) 



(13) 



(14) 



The above equations (|12I13I14|I show in each phase space 
a ffopf bifurcation of classical equilibria: we have a single 
point for the minimum energy Egs before the transition 
and an infinitely degenerate state represented in phase 
space by a circle. 

B. Chaotic regime: the equilibrium point below 
G+ /e = 1 bifurcates to two points ^23\ : 
i) trivial solution, stable for G+/e < 1 (and therefore 
G_/e < 1 ): 



91 = Pi = 92 = P2 = 0, (the origin) 



(15) 



ii) non-trivial solution I, stable for G+/e > 1 and G_/e < 
1 corresponding to the pitchforii bifurcation in each phase 
space: 

/2J(G2 -ew)V^' 
gi = 92 = 0, pi = ± ^ ^ j 

'J(G4 -£2^2) X 1/2 



P2 = T 



;2Gi 



(16) 



iii) non-trivial solution II, stable for G+/e > 1 and 
G_/e > 1: 



Pi = P2 = 0, gi = ± I ^ 



1/2 



92 = T 



/ J(Gi -eV 



2, ,2^\ 1/2 



(17) 



IV. ATOMIC WIGNER FUNCTIONS 

The maxima of the atomic Wigner distributions of the 
quantum ground state displays amazing similarity with 
the classical bifurcations: 

1) for A < Ac the atomic Wigner function (AWF) 
in both cases are a Gaussian-like state centered at the 
classical equilibrium position (the origin - not shown). 
2) for A > Ac, we separate the two cases: 
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Figura 2: (a) contour plots at the height of 50% of the In- 
tegrable Atomic Wigner functions of the ground state for: 
ijj = e = 1, \ = G/e — 1.5, G' — and scaled variables 
{qi/\/4J,pi/\/4J). Dot-dashed line - J = 4.5, dashed line - 
J = 10.5. The black continuous line - maximal region of the 
AWF, and the gray continuous line corresponds to the pro- 
jection of the minimum energy region of classical fixed points 
as given by Eg. 11311 . (b) AWF showing the negative part for 
J = 10.5. 

(a) Hopf bifurcation: for the integrable case, there is 
only one non-trivial regime (A = G/e > Ac), and we 
show in Figl^L a contour plot of the AWF on the scaled 
atomic phase space (qi / V4J , pi / \/4J) . These figures 
correspond to the ground state of the integrable system 
for J = 4.5 and 10.5 respectively at the same value of 
A = G/e — 1.5 > Ac- Note that the projection of the 
region covered by the maxima is an annular region for 
both cases, and is shown in black continuous line the 
circle corresponding to the maxima of the AWF. In gray 
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continuous line we show the region corresponding to the 
bifurcated classical stable fixed points. This shows that 
although qualitatively the mean field predicts correctly 
the change in the GS, the result does not give the 
position of the radius correctly. It is well known that 
mean field type calculations predict very well the energy 
(as shown for this model in '2^), but the states are 
not so well described. In Fig|2l3 a 3-dimensional plot is 
shown, with its negative part clearly visualizable for the 
case J = 10.5. 

(b) pitchfork bifurcation: for the non-integrable 
case, one regime is shown (A+ = G+/e > Ac and 
A- — G-/e < Ac) the non-trivial solution I is shown 
in FigEl Again, the maxima of the AWF clearly splits 
into a pair of peaks, connecting to the pitchfork type 
bifurcation shown in the classical case, and discussed in 

In order to estimate the number of basis states partici- 



0.8 p 

0.6 - 

0.4 - 

0.2 - 




5 



W ")ig4a" 
0.001 r 

0.0008 - . 
0.0006 - 
0.0004 - 
0.0002 - 




Figura 3: (a) Non-Integrable Atomic Wigner function of the 
ground state after the phase transition and the classical tra- 
jectories for G/e = G'/e = 0.75, J = 4.5 and E = -5.5, 
projected onto the atomic phase space with scaled variables 
(gi/\/4J,pi/\/4J); (b) another view of the same AWF show- 
ing the negative part. 



pating in the entanglement for different values of iV, we 
calculated the phase space area (in units of Ti) at the 
height of 50% of the AWF shown in FigHi for iV = 10 
and iV = 21. It becomes clear that it increases with N: 

A(J=4.5)^ o r-n J A(J=10.5)^ r or, 

thus a2i > aiQ. 

As for the non-integrable regime, Fig|5^ shows a very 
distinct AWF with two bumps centered at the bifurcated 
classical equilibrium points. We also plot the classical 
trajectories projected onto the atomic phase space for an 
energy slightly above the value of the GS. Note again the 
close analogy: the trajectories are concentrated in two 
distinct regions whose center correspond to the center of 
the two peaks in the quantum AWF. Fig|3j3 shows quan- 
tum interference effects, in the negative part of the AWF, 
but it is very small, not allowing to confirm the suspicion 
of ref. 28] that this GS is a Shrodinger-cat like state. It 
is also worth noting that there are quantitative as well 
as qualitative differences in what concerns the negative 
parts of the Wigner functions. In the integrable situation, 
its negative part is distributed in narrower region of phase 
space as compared with the non-integrable one, where 
negative contributions seem to be spread over a larger 
area. However, from the quantitative point of view, the 
integrable AWF exhibits a deeper negative contribution 



V. CONCLUSIONS 

In summary we have shown that the phase transition 
in Dicke's model can be qualitatively different in the in- 
tegrable regime and in the non-integrable one, the last 
case already discussed by other authors, although using 
the concurrence as a measure of entanglement. The lin- 
ear entropy produces the same qualitative results. 

The relationship between entanglement in the ground 
state at transition point, and the phase space bifurcations 
in the classical limit has already been noted in the liter- 
ature. We have found here an example of how the type 
of the quantum phase transition may qualitatively alter 
the corresponding classical instability. In the case of the 
Dicke model, the integrable and non-integrable situations 
are shown to be qualitatively very different. In quantum 
terms the GS exhibits a first order QPT reflected in the 
classical limit as a Hopf bifurcation which is very differ- 
ent from the already discussed non-integrable situation 
(second order QPT and pitchfork bifurcation). The close 
quantum and classical analogy is also strongly reflected 
in the corresponding GS atomic Wigner functions. 
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